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T-H ■ 

pj . An approximate partition functional is derived for the infinite-dimensional 

^ I Hubbard model. This functional naturally includes the exact solution of the 

Falicov-Kimball model as a special case, and is exact in the uncorrelated and 
atomic limits. It explicitly keeps spin-symmetry. For the case of the Lorentzian 
density of states, we find that the Luttinger theorem is satisfied at zero tem- 
perature. The susceptibility crosses over smoothly from that expected for an 
f^ ■ uncorrelated state with antiferromagnetic fluctuations at high temperature to a 

correlated state at low temperature via a Kondo-type anomaly at a characteristic 
temperature T*. We attribute this anomaly to the appearance of the Hubbard 
Q ■ pseudo-gap. The specific heat also shows a peak near T*. The resistivity goes 

^" , to zero at zero temperature, in contrast to other approximations, rises sharply 

around T* and has a rough linear temperature dependence above T*. 
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1 Introduction 

The single-band Hubbard model (HM) and the periodic Anderson model have attracted 
a great deal of interest since the heavy fermion and high temperature superconductors 
were discovered. Although there has been intensive theoretical work on both models, 
the knowledge we have for them is still limited, with an exact solution known only 
for the one-dimensional HM [|l|. An important advance was made by Metzner and 
Vollhardt , who observed that any diagrammatic treatment of fermion lattice mod- 
els simplifies in infinite dimensions. The proper self-energy becomes site-diagonal and 
momentum conservation is thus irrelevant [Q . Despite this simplification, the dynamic 
properties remain non-trivial. This discovery has stimulated much attention recently 
[^. Nearly all the phenomena associated with strong electron correlations have been 
found in the infinite dimensional models: the Mott transition |^, antiferromagnetism 
P, Fermi liquid or non- Fermi liquid behavior [0, and even the possibility of supercon- 
ductivity [§. The infinite dimensional models appear to be a good starting point for 
studying electron correlations. 

Brandt and Mielsch [Q first showed that the infinite dimensional {D°°) HM can be 
mapped onto an atomic problem in the presence of two time-dependent external fields. 
They were able to solve exactly the D°° Falicov-Kimball model, which corresponds 
to one time-dependent external field. A complete exact solution for the case of two 
external fields has not been possible. An equivalent mapping to the Anderson impurity 
model with a self-consistent condition has also been found |10, P). This has been used 



to calculate properties of the D°° HM [^ ^ and periodic Anderson model [|TT| by 
the quantum Monte Carlo (QMC) method. Although the QMC method is essentially 
exact, it can not treat the whole parameter region, in particular the cases of large U 
and low temperatures. Analytical continuation using the maximum entropy method 
to calculate some dynamic quantities is also problematic. A controlled approximate 
solution is clearly desirable for the whole parameter region. 

Many approximate schemes have been developed for the strongly correlated elec- 
tron systems, such as the alloy analogy approximation (plus coherent potential approx- 



imation), equation-of-motion decoupling |T2[, slave boson [|T3| and slave fermion [Q 



approaches, etc.. Unfortunately, they are usually not controlled approximations and 
do not become exact in physically-relevant limits. 

In the present paper we report on a better approximation scheme for the infinite 



dimensional models. The alloy analogy approximation |T2| (as well as the Hubbard-III 
approximation ||l5l and higher-order equation-of-motion decoupling ||T^) do rather well 
for high temperatures, and are exact in the free electron and atomic limits. But they 
violate the Luttinger theorem at T = 0, and do not include magnetic order probably 
as a result of breaking spin symmetry. The slave-boson and slave-fermion mean-field 
approximations usually give us contradictory results, and do not become exact even 
in the limit of free electrons. Our approximation scheme retains the merits of the 
previous approximations especially at high temperature and respects the spin symmetry 
explicitly. It satisfies the Luttinger theorem at zero temperature thereby remedying 
the most serious defect of approximations such as the alloy analogy approximation. 
Furthermore, the approximate solution is obtained in an expansion scheme so that it 



can be improved systematically by including higher-order effects. 

Janis and VoUhardt Jl^ have already suggested an approximate solution to the 



D°° Hubbard model. Their ansatz for the thermodynamic potential is the linear com- 
bination of the potentials for two Falicov-Kimaball models obtained by assigning one 
potential to one spin and retaining the other spin as the impurity. Although this ansatz 
is formally written to be spin-symmetric, the ingredients they used, i.e. the solution of 
the Falicov-Kimball model, break the symmetry intrinsicalUy. It does not satisfy the 
Luttinger theorem at zero temperature. 

We use the Brandt-Mielsch mapping to find an approximate solution to the 
partition functional of the D°° HM (including the periodic Anderson model) [|18] . From 



the partition functional one can obtain the self-energy and all the thermodynamic 
properties of the system as well as some correlation functions. The functional we 
obtain naturally contains the exact solution of the D°° Falicov-Kimball model as a 
special case, and is exact in the uncorrelated and atomic limits. The spin-symmetry is 
retained explicitly. 

To demonstrate the approximation scheme, we calculate the thermodynamic prop- 
erties of the D°° HM for a Lorentzian density of states (DOS). For this simple DOS, 
one can obtain some analytic expressions for physical quantities. The main physical 
properties of the D°° HM (except the Mott transition) are expected to remain for this 
DOS [^]. We find that at zero temperature the system is a Fermi liquid (except at 



half-filling), which is consistent with the exact result of Georges and Kotliar |T0[. The 
Luttinger theorem is satisfied. As the temperature is increased, the system crosses 
over continuously from a correlated state to an uncorrelated state with a Curie- Weiss 
susceptibility via a Kondo-like anomaly. This is in contrast to other approximations in 
which either the local moment behavior observed at high temperatures or the strongly 
correlated behavior observed at low temperatures is not properly reproduced |[l^ |T^. 



The Kondo-like anomaly here actually results from the Hubbard pesudo-gap. The spe- 
cific heat shows a peak at the corresponding (Kondo) characteristic temperature T*. 
The resistivity goes to zero at T = 0. This is completely different from the result of the 
alloy analogy approximation and other schemes, which break the Luttinger theorem at 
zero temperature so that the resistivity always goes up as temperature decreases. 

The paper is organized as follows. We briefly outline the main features of the 
D°° HM and the Brandt-Mielsch mapping in Section II. We derive the approximate 
partition functional by an expansion in Section HI. Some physical quantities, such as 
the susceptibility, the specific heat and the resistivity are then calculated using the 
functional in Section IV. Finally, we discuss our results in Section V. 

2 The infinite-dimensional Hubbard model 
and equivalent atomic problem 



2.1 The D^ Hubbard model 

The Hubbard Hamiltonian for large dimensions, D, is P| 

where Oj o-, a]^^ and n^ o- = al^^Ui^^ are usual fermion operators; fj, is the chemical poten- 
tial, E^ takes account of the spin-dependence of the atomic levels. In the presence of 
the external field h, E„ = Eq — aixsh {fJ^B the Bohr magneton). The hopping term is 
scaled so as to make the kinetic energy and the interaction term competitive. 

Because of the scaling by Xj \fD for the kinetic term in (|T]), the self-energy becomes 
site-diagonal [0, |^, i.e.. its Fourier transform S(iu;„,k) = S(iu;„) is a k-independent 
quantity. For a homogeneous system [^, the Dyson equation for the full momentum- 
dependent Green's function of (|I|) is given by 

G^{iun, k) = [iun - {E^ - /i) - S^(iu;„) - e{k)Y^. (2) 

Integrating G„{iuJni k) over k gives the following local Green's function 



po{e)de 



oo iUn - {E„ - /i) - T^aii^n) " e' 



(3) 



Here po(e) is the bare density of states (DOS), which is of Gaussian type 0] for the 
nearest neighbor hopping in (1): 



Po(e) = I] <^(e - Cfc) = ^expf-e^] 



The self-energy can be expressed as a functional of the local Green's functional G^°'^{iuJn)i 
as shown using the perturbation expansion in |pl| : 



S(z^„) = Y.[G':\iUn)]. (4) 

Once the self-energy functional ( ^) is known, one can determine the local Green's 
function from (|]), and, consequently, the self-energy. Eq. (|^) therefore provides a 
self-consistent equation for determining the Green's function of the system. 

2.2 The Brandt-Mielsch mapping 

The central problem is to find the self-energy, S, as a functional of the local Green's 
function G^°^{iujn)- Brandt and Mielsch mapped the 0°^ HM onto an atomic prob- 
lem with two external fields to find this functional. 

For an atomic (single site) problem, Hh in (1) reduces to the atomic Hamiltonian 
H = Ho + V, where 

Hq = ^(E^ - n)ala^, 
V = Un^n^. (5) 



Of course, the Hamiltonian H can easily be solved. However, when the two external 
fields are introduced, the problem is never trivial. 
A generating partition functional is defined as 



Z = Tr{e-^^^}, (6) 



where 



r/3 ri^ 



S = T.exp{- f dr f dr'Y.^''{r,r')a^\r)a^{r')} (7) 

JO JO „ 



with a^ (r) = e^'^aj.e"^'^ and \"(r,r') is an external field. The Green's function is 
defined by functional differentiation: Ga{T, r') = —6\nZ/6X''{T', r). Using Matsubara's 
formalism, we write 

^'^(^'O = ^T.^M-^^n{r-r')]X:, (8) 

P n 

where Un = {2n + l)7r//3. The Fourier transform of the Green's function, Gf^{iuJn), is 
then given by 

G.ituJn) = -d\nZ/d\l (9) 

The self-energy is defined via 

Y.„{iun) = G^^{iujn) -Ga{iuJnr, (10) 

where 

Gl{^u;n) = [^cOn-iE,-^^)-XT'■ (H) 

One may ask what the relation is between the local Green's function G^°'^{iujn) in (|^) 
and the Green's function G^iiUn) defined in (|). When one expands the self-energy of 
(|TUp in a series of skeleton diagrams in terms of G(j{iuJn), i-e., expands in the formalism 
of Baym and Kadanoff [^ , one finds that there is a one-to-one correspondance between 
these skeleton diagrams and those of the self-energy of the original Hubbard model in 
(1) in terms oiG^^'^{iUn)- The self-energy ([T0|), when written as a functional of (jo-(ico'„), 
therefore has the same form as the self-energy functional of the Hubbard model (1). 
Since this functional is what is inserted in (^, it is clear that G(j{iuJn) and G^^'^{iuJn) 
are the same once Eq. (H) has been solved self-consistently. 

In a real calculation of the local Green's function from the self-consistent equation 
@, one does not need to find an explicit form of the self-energy as a function oiGa{iuJn)- 
It is in general easier to determine from (^) the A^ (or equivalently the Green's function 
G^{iu!n) of (|llD ), and then to use the X" to fix the self-energy in ([T0|) and G^iiUn) in 
(^, which is equal to G^°'^{iUn). Below we will not differentiate between G^iiUn) and 
G^^^{iuOn) on the understanding that the external fields are determined self-consistently. 

A complete solution for the partition functional Z in (^ has not been possible. 
The exact solution is known for the Falicov-Kimball model 0, for which only one of 
external fields is needed p3[ as only one spin species is allowed to move in this model. 



In the next section, we derive an approximate solution for the partition functional Z 
of the D^ HM. 



3 Approximate partition functional 

3.1 Expansion 

Since Hq and V commute for the atomic Hamiltonian iJ in (|^), one can write 



oo 






e 



-/3^ori I l^-PU 



[l + (e-^'^-lKnJ. (12) 



The partition functional of (P) can therefore be expressed 

Z = e-^^n(5)o + {e-^'' - l)K^i5)o], (13) 

where VLq is the thermodynamic potential for Hq^ and (^)o = Ti{e^'^^°~^°^6). The im- 
portant point to note from (|13|) is that we need only to calculate expectation values 
with respect to Hq. However, the interaction still appears in the Heisenberg operators 
a^ (r) and a^ir) which define S* in (|^. 

We represent the operators a^ (r) and ^^(t) in terms of the Heisenberg operators 
for Hq in a similar way as in (|l^). We write 

af(r) = at(r)[l + (e^^-l)n_.(r)], 
a'^ir) = a.(r)[l + (e-^--l)n_.(r)], (14) 

here alir) = e^o^ate"-^o^. 



We first look at the expectation value {S)o in (|T^), which is given by 

f/3 rf^ 



{S)o = {Tre{exp[~ f dr j dT%{r,r')] 

Jo Jo 

exp[- [^ de [^ de%{e,e')]})o, (is) 

Jo Jo 



where 



; ' 7 5 



V^ir,r') = A^{r,T')B^{r,r' 

A^iT,r') = A'^(r,r')4(r)a.(r'), (16) 

B^ir,r') = [l + (e^--l)n_.(r)][l + (e-^-'-l)n_.(r')]. 

Expanding the exponentials in (llSD gives 



{S)0 = T.^ fdTfd9XmniT,9), (17) 

^ mini J J 

where 

XUr, 0) = (T.,{Af (r)5f (r)4(0)fif (^)})o, (18) 

with A^{t) = A^{T^,T[)...A^{Tm,T'J, and 5™(r) = 5,(ri, r()...5,(r^, r;;), and the 
integration in (|l^ is for all r's and ^'s. ^™(t) and B^^ij) contain all the operators 
involving spin-o" electrons in ([18D , and so using Wick's theorem we can write 



X„„(r, 9) = {Tre{Ar{T)Bl{9)}UTre{Al{9)B^{T)])o. (19) 



The quantity X^n may be calculated for any finite order m and n using Wick's theorem. 
However, an analytic calculation up to infinite order including all Feynmann diagrams 
as required in (|l^) is not possible. We therefore approximate to find the X^n- 

3.2 Approximation 

When evaluating the Xmn, we will consider only contractions within v4o-(t) and B^aiO) 
and ignore those between Acj{t) and B^cr{0)- Note that the effects of the hopping terms 
in the original Hamiltonian are incoporated via the external fields, A'^(r, r'), included 
in A(j{t) [see Eq. ([16|)]. Roughly speaking, the ^o-(r) account for the "dynamics" of 
the spin a electrons, while the B_^{6) monitor the value of rig. at the times, 6', at which 
—0" occupation changes. Restricting contractions to within B_„{6) sets up a mean '—a' 
field for the motion of the a electrons and vice versa. As we will see, (TqB^^IO)) gives 
rise to two terms which describe propagation in the upper and lower Hubbard bands. 
Our approximation is in the spirit of the Gutzwiller approximation for Gutzwiller 



wave function p4|. Gutzwiller's original idea was to treat the down-spin electrons as 
static when considering the motion of the up-spin eletrons, and vice versa. This ap- 
proximation does not break the spin-symmetry and keeps the translational invariance, 
so that Fermi-liquid behavior can exist at zero temperature. The ground state energy 
per site for the Hubbard model like Eq. (1) can be written as 

E = —Y.(l-''knka + Ud, (20) 

fc,(T 

where rika is the fermion distribution function and d is the double occupancy per site. 
The effects of one spin species on the other are incoporated into the renormalization 
factor, Qa. Eq. (pOD has been rederived statistically |2^ and the Gutzwiller approxima- 



tion has been shown to be exact in the limit of infinite dimensions [^ . 

Only considering the contractions within A„{t) and B_„{9) is equivalent to decou- 
pling A„{t) and 5_.(0). Writing (T.A-(r))o {TeBl„{e))^ gives 

n 

{TeBlM)o = {\l-n,]+^MUY.(e,-e[)}n„),. (21) 

Accounting for the first term of ([2l|), when multiplied by {T^A^{r))Q., is straightfor- 
ward. Summing over m for (T^A|"(r))o yields 

oo -| 
Wo = T.—^TrA^{T)), 



m=0"^i 



r/3 rP 



= (T,exp{- / dr [ dT'X'{r,r')aJ{T)a,{T')})o 
Jo Jo 

= exp{^ln[l-A;;^°(z^„)]}, (22) 

n 

with 



The second term provides a time-dependent "potential" [the exponential in (|2ll) ] 
which affects the dynamics of A"^(^) of the —a electrons. This term essentially changes 
the external field \_„{e,e') to A_^(^, ^')exp{t/(^ - 6')}. Summing over n for A%{e), 
one needs to calculate an expectation value of the type 

(S_^)o = (T,exp{- f de f dO'X-^'ie, ^')e^('-'')al,(0)a_.(0')})o. (24) 

Jo Jo 

These can not be dealt with as in (^) because of the non-periodicity of A^'^(^, 9')e'^^^~^ \ 
However, introducing the Hamiltonian 

and defining new Heisenberg operators with respect to H-a'- ctL^{9) = e^^a_^{6)^ we 
can easily calculate the expectation value 

Cs-.)h ^ = {exp{- Tde rde'x-'^ie,eraLm^Ao')})H ^, 

^ Jo Jo 

= exp{5:in[l-Ar£,(za;„)]}. (25) 

n 

Here 

~9-A^^n) = -. 7tA T^, (26) 

and {0)jj_ = Tre'^^^~''~^'^^9. It is clear that (^51) describes motion of electrons in the 
upper Hubbard band, while (0) comes from motion in the lower Hubbard band. This 
is reminiscent of the alloy analogy approximation (AAA), where the lower and upper 
Hubbard bands are introduced explicitly. However, as we shall see, our approximation 
avoids the most serious fault of the AAA and does not violate the Luttinger theroem 
at T = 0, which allows the metallic behavior. 

The expectation values of the operators B^^{6) always lead to a contribution pro- 
portional to ((1 — rio-)) and a time-dependent term {exp{U Yl^ii^i — ^Dl'^o-)- We 
then incorporate the exponential in the latter term directly into an operator A^^ 
= A"^^{9)exp{U J^^ii^i ~ ^'i)}- Ou^ decoupling procedure for this term is equivalent 
to writing for any function F[n-a{0)] 

T,{e^(^°-^°)i™,(0)F[n_.(0)]} 
= T4e^(^--~^-'')i™,(^)e^^"--F[n_,(0)]}e'^(^«-^--) 
^ {TeAZmH.Ae^^''-F[n.^{9)])s_y^'''^'''-\ (27) 

We can see that the upper Hubbard band is introduced naturally at the decoupling 
stage. This procedure must be carried out for -B"(t) as well. 
The final result for the decoupling of the Xmn in (|l^) is 

Xmn{T,9) ^ {TrA^ir))o{l-n^)o{TeA^imo{l-n^)o 

HTrA]^\r))f,y^-^{l-n,))f,/^^^^-''^\TeA-{e))o{n,)o 
+ (r.A-(r))oK)o(T,ij(^))^^(e^^"Kl-^i)^,e^(^"-^^) 

+ Mr(^))i/T<^'''"'^t)H,(T.^r(^))H,(e'^^"^ni)H,e'^^^^^"'^^^(28) 



where Q^ is the thermodynamic potential of Her. Eq. (|28|) is the main approximation 
we make for the partition functionaL 

This approximation leads us to a form for {S)q 

{S)o = (5t)o[l - K)o] (5i)o[l - K)o] 

+ {S,)H/''{n,)H,{S,)f,/^{n,)ji/^^^''^-^^\ (29) 



where {Sa-)o and {Scr)H„ is given by (^) and (pSj), respectively. The expectation value 
of the term (ri|?7,|S')o in ([13D is then given by 

(ntn;5)o = (ST)^^e^^ (nT)H,(^i)H,e^^(n,)^^e^i:.(^-^^). (30) 

Using ([13|), (p9|) and (pOD, we obtain the partition functional 

Z = ( Z, Zt ) f l; ""'!<' " "'' "^^' " "'* ) ( f ) ■ (31) 

Here 

Z^ = exp{^ \n[iujn - {E^ - yu) - A^]}, 

n 

Z. = exp{^ln[z^„-(E.-/i)-f/-A^]}, 

n 

n^ ^ {n^), = [e^iEa-f^) + lf' (32) 

The partition functional Z in ( pi] ) has a very natural interpretation, n^j and Z^- can 
be regarded, respectively, as the "occupancy" and the "partition function" of the spin 
a particles in the lower Hubbard band, while "h^ and Z„ are the corresponding quan- 
tities in the upper Hubbard band. The terms with {1 — n^){l — ni) , n|(l — fij) (or 
(1 — n'\)ni ) and e^^h-^ni in the matrix account for holes, singly-occupied and doubly- 
occupied configurations, respectively. The factor e^^ is essential because it cancels the 
double counting of the interaction U involved in Z-^Z^. Comparing the form of the 
ground state energy E of the Gutzwiller approximation in (|2D| ) with Z in (PTf), we see 
that the "configuration" matrix in Z gives a function like the renormalized constant 
Qa in E. 

For the D°° Falicov-Kimball model, where one of the external fields, A^ say, equals 
zero, Eq. ( pi] ) reduces to 

ZpK = Z^ + Z^exp[-P{E^ - fi)], (33) 

which is just the result obtained by Brandt and Mielsch 0. Therefore, our partition 
functional naturally includes the exact solution of the Falicov-Kimball model as a 
special case. It is also easy to check that ( pT] ) gives the exact result for the atomic 
(t = 0) and ?7 = limits. 

9 



4 Thermodynamic properties 

4.1 The Green's function and physical quantities 

Using (^) and ([To|), one can obtain the self-energy functional. Substituting this func- 
tional into the right hand side of the self-consistent equation in (^, and noting from 
@ that the left hand side can be written 

G':\iuOn) = V^n - {E. - m) - A^ - S.(2cu„)]-\ (34) 

one can determine the external fields self-consistently. Consequently, the self-energy 
and the Green's function for the D°° HM are easily found using this approximation. 

If the external fields are fixed, it is convenient to use another form of the local 
Green's function derived directly from the functional derivative in (P), 



where 



Gf («^n) = Ga{^UJn)Po + QA^^n)Pl (35) 



G^{iuJn) = [iuJn - {E„ - ;u) - A;^] \ 

G„{%u:n) = [tcOn-iE^~^)-U-\X\ (36) 



and 



^0 = -[Z^{l-n^){l-ni)Zi + Z^n^{l-hi)Zi], 
Pi = |[^T(1-^T)%^T + V^%^i^i]- (37) 

Clearly, Pq + Pi = 1. Again, the Green's function G^°^{iuJn) looks similar to that given 



by the alloy analogy approximation |]12| . However, Pq and Pi have quite different forms. 



We will see later that this difference leads to important consequences. 
Using (|35|) , the density states of the system can be expressed as 

D^itUJ^) = V^itUJn)Po + 'Da{l^n)Pl (38) 

where 

vs^) = -^[g„{uj + to+)-g^{uj-to+)], 

Ztxi 

V^iiun) = -:—[ga{uJ + tO^)-g^{uj-tO+)]. (39) 

Zm 

The chemical potential should be determined self-consistently from the following equa- 
tion for the electron concentration 



n,{T) = / def{e)D^{e), (40) 

where /(e) is the fermion distribution function and D„{e) is given in (^). 

10 



From the Green's function (|35|) , one can calculate all thermodynamic quantities of 
the systems. In this paper we discuss the suceptibility, specfic heat and resistivity of 
the system. 

The static susceptibility can be obtained directly by differentiating the occupancy, 
no-, in pop with respect to Ea-. We set E„ = —a fish, with h the external magnetic 
field. The susceptibilty, xiT), is then 

X{T) U=o- -2/iB(^ - ^) \e,=e,=o (41) 

In order to obtain the specific heat, it is convenient (especially for numerical cal- 
culations) to start from the internal energy, E{T). This enables us to calculate the 
specific heat C{T) as the first derivative of E{T) with respect to T rather than the 
second derivative of the free energy. In the infinite-dimensional case, E{T) is given by 

m 

Em = - E / *po(.) / ^/(.)(. + .)lu.-^^^-^l-__ + 1,„, (42) 

where Po(e) is the bare density of states. 

The conductivity is expressed using Kubo's formula as a current-current correlation 
function. Usually this correlation function can not be calculated exactly since the ver- 
tex corrections are difficult to treat. However, in the infinite D limit, vertex corrections 



vanish [^, so that one is left with a simple bubble for the conductivity correlation 



function. The conductivity becomes p6| 



oo roo 



a(0) = 2/ depoie)! ^(-^M(^,^)^ (43) 



oo v/— oo 



where the spectral function A{e,uj) is defined as 

A{e,uj) = i[Grie,uj) -Ga{e,uj)] 

1 1 



u;-|-yU — Sr + zr/ — e u + fi — I^a — irj — e 



(44) 



with Sj, and S^ the retarded and advanced self-energy of the system, respectively. We 
have set h = e = a = 1 (a is the lattice constant) in (^31). 



4.2 The Lorentzian bare density of states 

To illustrate the scheme we use a Lorentzian type for the bare DOS []TU|, |T^: 

Po(e) = -(e^ + n. (45) 

IT 

Although this DOS has some undesirable properties as a result of having too much 
weight at large energies, it leads to some analytical results and some analytical expres- 
sions for physical quantities. Apart from the Mott transition, the physical properties 



are expected to be qualitatively correct for this DOS [jT9|. 
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For the Lorentzian DOS, Eq. (H) reduces to 

G':\iUn) = [luj - {E, - /i) - S,(z^„) + iTsgnUn]-\ (46) 

Comparing this with Eq. (0), one immediately obtains the external field 

A„ = -iVsgnUn- (47) 

Substituting (^71) into (^) and ([T0|) gives the Green's function G^°^{iuJn) and the self- 
energy Ti{iuJn) for a system with a Lorentzian DOS. A frequency summation then gives 
Z„ and Zfj in (|3T|) 

Z, = C'e^"-(^'^+^-^\ (48) 

where C is a constant, which is irrelevant for the thermodynamic properties, and the 
function 

/• — X /"OO 

u^{x) = dX f{z - X)poiz)dz. (49) 

Jo J-oo 

At zero temperature, Z^- in (|^) reduces 

Z. = Cexp{/3{^^^ + -[(/x-iJ;.-[/)arctan^^;^-^ln(l + ^~f^' )]}}, (50) 

while Zcr has the same form as Z^r but with E^r replaced by E^r + U. 

We find that for the electron concentration n {= n-^ + n|) less than one, the term 
with (1 — n^){l — nj.) in Z of (^) is dominant, and Pi vanishes. The chemical potential 
can easily be obtained from (HO) 



/i = — rtan— (1 — n) (for n < 1). (51) 

For n > 1, the term with e^^n-^n^ in Z is dominant (so that Pi = 1) and the chemical 
potential is given by 

/i = f/ + rtan-(n-l) (for n > 1). (52) 

For n = 1, the terms with n^{l — n^) [and (1 — ni)ni ] in Z are important and the 
chemical potential fi = U/2, as it should be. The chemical potential n is shown in Fig. 
1 as a function of the electron concentration n. 

From (0) we find that the imaginary part of the self-energy at cu = is equal to 
zero away from half-filling and the Luttinger theorem is satisfied. The system is in a 
Fermi-liquid state. This is consistent with the exact statement of Georges and Kotliar 
[ PHI . However, at half- filling the imaginary part of the self-energy at to" = is not equal 
to zero and there is a discontinuity at the chemical potential, as shown in Fig.l. At 
half-filling the system is therefore a Mott insulator for any finite U. This is similar 
to what happens in the one-dimensional case |1|, but is inconsistent with Monte Carlo 
calculations |^, which show that the Mott insulating state appears for all U >3. The 
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occurrence of a Mott insulating state at half-filling follows in our approximation as 
a result of the particle-hole symmetry and is not an artefact of using the Lorentzian 
density of states. The absence of a Mott transition at finite U in our approximation 
results from overcounting of the interaction effect for small f/, i.e. for the same reason 



as given by Hubbard for his Hubbard-II approximation p8 
The susceptibilty is given by 



2 r-Po r Pq 



which has no divergence for any n and finite [/, so that there is no ferromagnetic phase. 
The temperature dependence of the static susceptibility x and the specific heat 
C are shown in Fig.2. (The temperature T and the interaction U are in units of 
r) At high temperature, x shows a Curie- Weiss-like behavior (x~l/(T + 0) with 
^ > 0), indicating that the system behaves like a system of independent moments with 
antiferromagnetic fiuctuations. For n near half-filling or for large [/, x decreases as T 
drops below a characteristic temperature, T*, in a way reminiscent of the Kondo effect 
(see curves A and B). Near T*, the specific heat also has a peak. For the susceptibility 
there appears to be a second temperature Tj*, below which a correlated state is formed 
where the susceptibility is nearly constant, i.e., Pauli-like. At T^ there is no anomaly 
in the specific heat. For small n, the susceptibility continuously crosses over from 
the Curie- Weiss behavior to the Pauli-like behavior without any Kondo-like anomaly. 
Correspondingly the specific heat shows no peak. At very low temperature, the specific 
heat as a function of T is consistent with a power law, but with an exponent not exactly 
equal to one, i. e. the behavior is not exactly linear. The case of the Lorentzian density 
of states corresponds to the Anderson impurity model |TD[ with e^ = /i(T), where e^ 



is the d-level energy. The impurity model has been solved exactly for fixed e^ and 
does give C„ oc T. However, in our case the chemical potential /i is a function of 
temperature, which must be determined self-consistently. This has not yet been done. 
The variation of /x with T will change the impurity model result for C^. It is therefore 
not clear whether our result is actually correct or a result of the approximation. 

In order to understand the "Kondo anomaly" in this system, we calculate the den- 
sity of states, shown in Fig. 3 for the case of f/ = 2 and n = 0.8. There are two 
peaks above a certain temperature, which corresponds to the "Kondo" characteristic 
temperature T*. There is only one peak below this temperature. The Kondo anomaly 
seems to be related to the appearance of the Hubbard pseudo-gap. The self-consistent 
perturbation calculations also show one peak in the DOS of the system at zero tem- 
perature |21|, ^. However, using the perturbation scheme of Yosida and Yomada ||3CI|| , 



Georges and Kotliar |jTO| found that the DOS has two peaks near half-filling at zero 
temperature. The one-peak feature found at very low temperature obtained using the 
self-consistent pertubation calculation might indicate that the upper band does not 
contribute much to thermodynamic properties at low temperature and the Hubbard 
pseudo-gap behavior manifests itself only above a characteristic temperature. 

The temperature dependence of the resistivity p (= l/o") is shown in Fig. 4 (for U = 
2) and Fig. 5 (for U = 4) for various electron concentrations. The resistivity goes to zero 
when temperature decreases to zero as it should, given that our approximation satisfies 
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the Luttinger theroem at zero temperature. This is in contrast with the result of the 
alloy analogy approximation (shown shematically by the dashed curve in Fig. 4), which 
keeps rising as temperature decreases. The resistivity becomes roughly linear with T 
above the characteristic temperature, while it rises rapidly around the characteristic 
temperature. In order to see clearly how fast p changes with T, we show in Fig. 6 the 
first derivative of p with respective to T in the case of f/ = 2. There is a sharp peak 
at the characteristic temperature. 

5 Discussion and conclusion 

In this paper we have derived an approximate self-energy functional for the infinite- 
dimensional Hubbard model. It naturally includes the exact solution of the D°° Falicov- 
Kimball model as a special case and retains the spin symmetry. This finite-temperature 
theory successfully incorporates the high-temperature uncorrelated behavior and the 
strongly correlated behavior at low temperature in a unified way. Many approxima- 
tions used before, such as the alloy analogy approximation, the equation-of-motion 
decoupling, the extension of the Gutzwiller approach and the slave boson mean-field 
theory, can not reproduce both limiting behaviors correctly without introducing a spu- 
rious phase transition |T^ |T^. It seems to us that our approach is the first to give 
physically reasonable results for all temperatures. 

The functional derived here is very easy to use. It can also be directly applied to the 
periodic Anderson model. The only change is to replace the self-consistent equation in 
©by 

Gt = r '-'^ ^^, (54) 



iujn - {e- p) - Ef^{iuJn) - 



iun - (e - /i) 

where "/" represents localized electrons, and V is the mixing energy between conduc- 
tion and localized electrons. From (Q) one can determine the local Green's function 
and the self-energy of the localized electrons, and consequently, all Green's functions 
of the system [^ . 

In this paper we have treated only the case of a homogeneous system. For some 
symmetry-broken phases, such as the antiferromagnetic state, the self-consistent equa- 
tion changes |^, although the form of self-energy functional remains. Becuase our 



scheme treats spin-up and spin-down electrons on an equal footing (in contrast to the 
alloy analogy approximation), we expect that it can describe correctly the antiferro- 
magnetic instability. We will discuss this issue elsewhere |^ . 



There exist some weaknesses in the present form of the partition functional. The 
dynamic flucuations induced by the motion of one spin species to another are not prop- 
erly treated. Moreover, the Mott transition does not take place at a finite interaction 
U . Since our approximation is an expansion, we can include relevant higher order con- 
tributions which account for these flucuations. The Mott transition has been produced 
by Hubbard in his Hubbard-III paper [15|| by including the resonant correction in the 



equivalent coherent-potential approximation. It would be interesting to examine this 
transition in our formalism by including higher order fluctuation effects since our solu- 
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tion satisfies the Luttinger theorem in the inetahc side at zero temperature, while the 



Hubbard-III solution does not 133 
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FIGURE CAPTIONS 

Fig.l. The chemical potential /i as a function of the electron concentration n at zero 
temperature, /x is in unit of F. The gap at half- filling is present for all non-zero 
U 

Fig. 2. The temperature dependence of the susceptibility, x (solid curve), and the 
specific heat C (dashed curve) for the Lorentzian density of states for various n 
and U . Curves A &i a are for U = 1 and n = 0.9, while B h h^ C h c and D & 
d correspond to n = 0.8, 0.6 and 0.4, respectively, for the case of f/ = 2. T and 
U are in units of F, and x is in unit of /i^. 

Fig. 3. The density of states of the system, D{uj), for various temperatures in the case 
of f/ = 2 and n = 0.8. 

Fig. 4. The temperature dependance of the resistivity for U = 2 and various electron 
concentration n. The result of the alloy analogy approximation is schematically 
shown by the dashed curve. 

Fig. 5. The temperature dependance of the resistivity for U = 4 and various n. 

Fig. 6. The first derivative of the resistivity with respect to temperature for U = 2, 
corresponding to Fig. 4. 



